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1.  Introduction 


The  increased  interest  in  6-systems  during  the  recent  years  (see  [T6j,  and  references 
therein)  is  due  mainly  to  two  reasons:  (a)  6-systems  provide  superior  finite  wordlength 
properties  with  respect  to  roundoff  noise  propagation  [5]  and  coefficient  sensitivity  [l],  [5], 
(7),  as  compared  to  their  ^-system  counterparts,  and  (b)  the  6-operator  yields  the  differ¬ 
ential  operator  as  a  limiting  case  when  sampling  time  approaches  zero  enabling  a  unified 
treatment  of  both  continuous-  and  discrete-time  systems  [1]. 

With  such  advantages  in  mind,  development  of  2-D  and  multi-dimensional  (m-D)  6- 
system  models  must  clearly  be  undertaken.  Such  research  can,  for  example,  provide  m-D 
digital  filters  with  superior  roundoff  error  and  coefficient  sensitivity  performance  allowing 
their  implementation  to  be  carried  out  in  a  shorter  wordlength  environment.  This  is 
especially  crucial  in  real-time  applications,  such  as,  in  implementing  narrow  bandwidth 
filters  under  high  sampling  rates  (for  example,  in  current  wide  bandwidth  communication 
system  applications)  where  traditional  ^-operator  implementations  perform  poorly  [8]. 

In  applications  mentioned  above,  and  those  dealing  with  high-speed  processing  of  2-D 
and  m-D  data  (for  instance,  in  weather,  seismic,  gravitational  photographs,  video  images, 
systems  with  mutliple  sampling  rates,  etc.),  ensuring  stability  is  an  important  consideration 
(see  [9],  and  references  therein).  Given  the  characteristic  polynomial  of  a  6-system,  to 
determine  stability,  one  may  first  use  a  variable  transformation  that  yields  a  more  familiar 
stability  region,  for  instance,  the  unit  bi-circle.  Then,  an  existing  technique  (see  [9-10], 
and  references  therein)  may  be  applied.  However,  such  techniques  are  known  to  be  prone 
to  numerically  ill-conditioning  [1],  [6].  In  the  1-D  case,  direct  stability  checking  methods 
for  6-system  polynomials  are  in  [6]  (where  a  tabular  method  based  on  the  work  in  [11] 
is  given)  and  [12]  (where  a  Hermite-Bieler-like  Theorem  is  utilized).  Hence,  our  purpose 
here  is  to  develop  a  direct  easily  implementable  stability  checking  technique  applicable 
to  m-D  6-systems.  As  usual,  for  notational  simplicity,  we  concentrate  on  the  2-D  case,  the 
extension  to  the  m-D  case  being  quite  straight-forward. 

In  checking  stability  of  bivariate  characteristic  polynomials,  two  conditions  must  be 
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satisfied. 

(a)  Condition  I  involves  a  1-D  stability  check  of  a  polynomial  with  real- valued  coefficients. 
One  may  use  the  table  form  in  [6].  Alternately,  one  may  utilize  an  explicit  root  location 
scheme. 

(b)  Condition  II  involves  a  stability  check  of  a  polynomial  with  complex- valued  coefficients 
where  the  latter  are  dependent  on  a  parameter  taking  values  on  a  certain  circle  in  the 
complex  plane.  Explicit  root  location  schemes  are  now  ineffective,  and  the  value  of  tabular 
methods  becomes  apparent.  Note  that,  in  such  a  situation,  compared  to  Nyquist-like 
techniques  [13],  tabular  methods  are  known  to  provide  certain  numerical  advantages  as 
well  [14]. 

In  checking  condition  II  for  2-D  ^-systems,  an  effective  technique  involves  checking 
positive  definiteness  of  the  Hermitian  Schur-Cohn  matrix  [15].  This  lets  one  use  an  impor¬ 
tant  simplification  due  to  Siljak  [16].  The  tabular  form  in  [15]  takes  full  use  of  this  since  it 
provides  the  Schur-Cohn  minors  (that  is,  the  principal  minors  of  the  Hermitian  Schur-Cohn 
matrix)  directly  from  its  entries  [15],  [17].  A  similar  simplification  applicable  to  6-systems 
is  clearly  possible  if  condition  II  may  be  reduced  to  checking  positive  definiteness  of  a 
Hermitian  matrix. 

With  the  above  in  mind,  we  develop  the  following  in  this  paper:  (a)  Tabular  form 
for  stability  checking  of  6-system  characteristic  polynomials  possessing  complex-valued 
coefficients,  (b)  Analogs  of  Schur-Cohn  minors  and  a  corresponding  Hermitian  matrix 
applicable  for  such  systems,  and  (c)  a  direct  stability  checking  algorithm  for  2-D  6-system 
transfer  function  models. 

The  paper  is  organized  as  follows.  Section  2  introduces  the  notation  used  throughout 
and  a  brief  review  of  previous  results.  Section  3  develops  a  tabular  form  for  stability 
checking  of  6-systems  with  complex- valued  coefficients  and  some  important  relevant  results. 
Section  4  presents  quantities  that  may  be  regarded  as  the  analogs  of  Schur-Cohn  minors 
for  6-systems.  The  2-D  stability  checking  algorithm  in  Section  5  is  based  on  the  tabular 
forin  for  real-valued  coefficients  [6].  Since  only  little  extra  work  is  needed,  results  in  both 
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Sections  3  and  4  however  are  developed  for  the  more  general  complex-coefficient  case. 
Section  6  presents  an  example  to  validate  the  results.  Section  7  contains  the  conclusion 
and  some  final  remarks. 
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2.  Preliminaries 


2.1.  Notation 


3?,  5  Real  and  complex  number  fields. 


3?px<,1  3pX?  Set  of  matrices  of  size  pxq  over  9?  and  respectively. 


var{-} 

R*N»  ImH 
f] 

At,  A ,  A* 
&Hn,  3tu>]n 

»(w) 
"l  K1 


3?(wi,W2) 


Number  of  sign  changes  in  the  sequence  {•}  of  real 
numbers. 

Real  part  and  imaginary  part  of  [•)  G  5. 

Complex  conjugate  of  [•]  G  §. 

Transpose,  complex  conjugate,  and  complex  conjugate 
transpose  of  A  G  9?x,)  respectively. 

Set  of  univariate  polynomials  of  degree  n  (with  re¬ 
spect  to  the  indeterminate  w  G  9)  over  9?  and  S, 
respectively. 

Set  of  rational  univariate  polynomials  (that  is,  quo¬ 
tient  of  univariate  polynomials)  over  9?. 

Set  of  bivariate  polynomials  of  relative  degrees  ni 
and  n2  (with  respect  to  the  indeterminates  ui  G  3 
and  tu2  G  O',  respectively)  over  9 R. 

Set  of  rational  bivariate  polynomials  over  9?. 


z,  c  Indeterminates  of  q -  and  5-systems,  respectively. 


r  Real  positive  number,  usually  the  sampling  time. 


The  transformation  relationship  between  corresponding  q- and  5-systems  is 


5  = 


q  - 1 

T 


Z  -  1 
r 


H 


H(c)  |c_, 
G(*)U- c 

H(Cl ,  C2  )|c— >z 

G(zi,z-i)\t^c 


♦/-system  quantity  analogous  to  its  corresponding  5- 
system  quantity  [•];  for  example,  transfer  function 
of  a  given  discrete-time  system  is  either  H(c)  if  im¬ 
plemented  based  on  the  5-operator  or  H(z )  if  imple¬ 
mented  based  on  the  ^-operator. 

H{c)\c=(z  —  \)/t 

G(z))z=l  +  rc 

H{ci ,  c2 ) | Ci.  =( * _ i )/r<  ,= j  )2 
G(zi,Z2)\zi  =  l  +  TCi,i=l,2 


(2.1) 
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Stability  studies  of  1-D  and  2-D  q-  and  6-systems  involve  the  following  regions: 


U] 

{zG 

3 
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{(c,,c2)e3!: 

|c,  +  1/t[  = 

l/r. 

t  = 

1,2}. 

To  avoid  unnecessary  notational  complications,  the  sampling  time  in  both  horizontal  and 
vertical  directions  is  taken  to  be  equal  to  r  >0. 

To  emphasize  the  degree  of  F(w)  =  Y^k=oa e  3[u>]n.  we  sometimes  denote  it 
as  F(w)n  as  well. 

F( to)  Conjugate  polynomial  of  F(w),  that  is, 

F*(z)  Reciprocal  polynomial  of  F(z),  that  is,  znF(l/z) 

F*(c)  Reciprocal  polynomial  of  F(c),  that  is,  (1  4-  rc)n 

*(n*0 

A  g-system  polynomial  is  q-symmeiric  if  F(z)  =  F^z).  A  6-system  polynomial  is  6- 
symmetric  if  F(c)  =  F*(c). 


Tabular  forms  of  stability  checking  of  a  polynomial  in  3[u>]n  typically  employ  a  sequence 
of  polynomials  each  of  descending  order.  The  first  row  of  such  a  tabular  form  is  denoted 


as  row  #n,  the  second  row  is  row  #n  —  1,  and  so  on. 

JT,  MJT  Jury  table  [18],  modified  Jury  table  [15],  [17], 


real-^-BT 

complex-(?-BT 

real-6-BT 

complex-6-BT 


Bistritz  table  for  <7-system  polynomials  with  real¬ 
valued  coefficients  [11]. 

Bistritz  table  for  <7-system  polynomials  with  complex¬ 
valued  coefficients  [19]. 

Table  form  for  6-system  polynomials  with  real-valued 
coefficients  [G] . 

Table  form  for  6-system  polynomials  with  complex¬ 
valued  coefficients  (this  paper). 
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A  ^-system  polynomial  with  all  its  roots  in  Uq  (for  the  1-D  case)  or  Uq  (for  the  2-D  case) 
is  said  to  be  stable.  The  corresponding  regions  for  a  (5-system  polynomial  are  U 4  (for  the 
1-D  case)  or  iXf  (for  the  2-D  case),  respectively. 


2.2.  Review  of  complez-q-BT 

The  complex-<5-BT  introduced  in  Section  3  is  based  on  the  complex-g-BT,  and  hence,  we 
briefly  review  it  now.  For  more  details,  see  {10].  Let  the  characteristic  polynomial  of  a 
(/-system  be 


F(z)  =  ^4">2*  €  $»(z|„  with  F(l)  €  &  and  F(l)  £  0. 


(2.2) 


k=0 


The  complex-g-BT  is  formed  using  the  symmetric  polynomial  sequence  (T(z), }’l=0 
where  [19] 


(  F(z)„  +  F‘(z)„, 


T(z)i  =  / 


for  t  =  n; 

for  t  =  n  —  1;  (2.3) 

[  (Si^  +  ji±2£gl£!i±i  -  rMit»  ,  fori=n-2,n-3,...,0, 


-  F»(^)„ 

Z  —  1 


where 


t 


X  -  _  to _  ■  _  _  n  _  q  n 

i+2 ~  f(o)i+I  “  #+■> ’  ’  ’ 


(2.4) 


As  in  [11]  and  [19],  equating  similar  powers  on  either  side,  we  may  also  get  the  following 
determinental  rule:  For  k  =  0, 1, . . . ,  z,  and  i  =  n  —  2,  n  —  3, . . . ,  0, 


-  _L_ 

l0 


p.i+2) 

c0 

i<H-2) 

H+1 

1 

1 

i<*+2) 

*•1+2 

A*+2) 
lk+ 1 

l0 

l/t+i 

+  *«+0 
U'+i 

U+i 

jU+1) 

lk 

^lk+ 1  • 


(2.5) 


Remark.  The  computational  advantage  of  BT  is  due  to  T(z),  being  ^-symmetric.  This 
(»’) 

'■k  *' 


implies  =  f,_t  ,  k  =  0, 1, . . .  ,  i,  and  hence,  it  is  necessary  to  evaluate  only  half  the 


coefficients  of  each  row. 

Using  (12-13),  (16),  and  Theorem  6  of  [19],  we  get 


ft 
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THEOREM  2.1.  [19]  The  polynomial  F(z)  €  9‘[z]n  is  ^-stable  iff 

I.  ^  0,  *  =  n  —  1,  n  —  2, . . . ,  0,  and 

II.  i/n  =  var{f(l)u,f(l)n_1,...,f(l)o}=0. 

2.S.  Some  results  on  2-D  stability 

Consider  the  2-D  ^-system  transfer  function 

H(zi ,  z2)  =  jr*1-’  2~;  €  ,  z7 )  (2.6) 

F(zuz2) 

where  F(zi,z2)  £  $[zi]nt[£2]ns  and  F(zi,z2)  £  ft[zi]m[z2]ni.  The  2-D  z-transform  is 
taken  using  positive  powers  of  Z{.  For  a  comprehensive  discussion  regarding  stability  of 
such  systems,  see  [9-10],  and  references  therein.  Hence,  for  reasons  of  brevity,  only  some 
analog  results  applicable  to  2-D  6-systems  are  provided.  It  is  only  necessary  to  observe 
that  the  corresponding  6-system  /7(cj ,  c2)  satisfies 

ff(c,,c2)  =  =  #(zi,*2)|,-.c  €  *(ci,c2)  (2.7) 

r(ci,c2) 

where  E(ci,c2)  £  3?[ci]ni [c2]„2  and  F(ci,c2)  £  8?[ci]ni [c2]nj-  In  the  remainder  of  this  pa¬ 
per,  we  will  only  be  dealing  with  transfer  functions  H(c\ ,  c2)  that  are  devoid  of  nonessential 

singularities  of  the  second  kind  on  Tf  and  the  pair  E{c\,c2)  and  F(c!,c2)  is  taken  to  be 

_ 2 

coprime.  If  the  2-D  polynomial  F(ci,  c2)  ^  0,  V(ci,c2)  £  Ug,  it  is  said  to  be  6 -stable.  After 
using  (2.1),  the  following  result  follows  directly  from  [20]: 

THEOREM  2.2.  The  2-D  6-system  in  (2.7)  is  6-stable  iff 

I.  F(cj,  -1/r)  ^  0,  Vc!  £  Us,  and 

II.  F(ci,c2)  7^  0,  Vci  £  Ts,  Vc2  £  Us- 

The  following  result,  which  allows  one  to  use  the  real-6-BT,  is  directly  from  [21-22] 
after  using  (2.1): 

THEOREM  2.3.  The  2-D  6-system  in  (2.7)  is  6-stable  iff 
I.  F(ci ,  -1/r)  ^  0,  Vci  £  Us,  and 
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II.  G(x,c2)  ^  0,  Vx  €  (— 2/r, 0],  Vc2  e  U6. 
Here  G(x,c2)  =  F(cx ,  c2)F(c, ,  c2)  ei6T< 

*=*<«I  +  C| 

2.^.  SchuT-Cohn  minors 


In  stability  checking  of  2-D  ^-systems,  the  following  result  is  important: 

THEOREM  2.4.  [15],  [23-24]  The  polynomial  F(z)  6  $[z]n  is  stable  iff  A,  >  0,  i  = 
1,2, ...  ,n,  where  A^  is  the  principal  minor  of  the  Hermitian  Schur-Cohn  matrix  r  =  f  *  = 
{7i>}  €  SnXn  defined  as 

« 

7 ij  =  ^  '(an— i+kan— j+k  ai—kaj—k)i  f°r  1  —  j- 

k=l 

Stability  checking  of  2-D  9-systems  then  involves  positivity  checking  of  all  Schur-Cohn 
minors  A,(z),  Vi  =  1,2, ...,n,  V|z|  =  1.  A  necessary  and  sufficient  condition  for  this  is 
positivity  of  Ai(l),  Vt  =  1,2, ...,n,  and  An(z),  V|z|  =  1.  This  is  the  simplification  due 
to  [16]  that  has  been  effectively  utilized  in  applying  the  MJT  [15].  The  advantage  of  the 
latter  is  that  its  entries  yield  the  Schur-Cohn  minors  directly.  The  fact  that  complex-g-BT’s 
entries  also  yield  the  Schur-Cohn  minors  was  only  recently  shown. 

THEOREM  2.5.  [i0],  [25]  The  Schur-Cohn  minors  of  F(z)  are  the  principal  minors  of 
the  (n  x  n)  tridiagonal  Hermitian  matrix 


A  = 


illi-V'ii’-5'! 
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3.  Complex-6-BT 


With  no  loss  of  generality,  consider  the  ^-system  characteristic  polynomial 


where 


Oq  (=  and  ag  >  0. 


(3.1) 


(3.2) 


We  now  construct  the  complex-£-BT  with  the  use  of  the  ^-symmetric  polynomial  se¬ 
quence  {T(c)i}"_0  where 

(F(c)n  +  F«(c)n,  t  =  n; 


T(c)i  = 


[cU-t’lc 


",  *  =  n  -  1; 


Here 


(fr+2  +  6j+ 2(1  +  rc))r(c),-n  -  T(c)i+2 

V  1  +  TC  ’  1  -  n  ^ 

_  T(-l/r)i+2  n 

5,+2_  T(-l/r)i+1’  ’  ’ 


The  normal  conditions  required  to  complete  the  sequence  are 

T(-l/r)^0,  t= 


(3.5) 


Remarks. 


1.  To  determine  ^-stability  of  F(c),  one  may  of  course  first  obtain  F(z )  =  F(c)|c-_*  and 
then  determine  its  ^-stability  by  applying  familiar  stability  checking  algorithms  (e.g.,  BT 
or  MJT).  The  possible  shortcomings  of  such  a  scheme  are  outlined  in  [1]  and  [6].  The 
purpose  here  is  to  obtain  a  direct  check  for  ^-stability. 

2.  We  follow  the  work  in  (6)  and  [19],  and  hence,  for  brevity,  all  details  are  omitted. 

3.  The  conditions  T(— 1/r),  =  0,  for  some  i  =  1,2, . . .  n  —  1,  imply  certain  singular  condi¬ 
tions  on  the  root  distribution  of  F(c)  [11],  [19].  The  equivalent  singular  conditions  for  the 
real-<5-BT  is  in  [6]. 

4.  Using  ^-symmetry,  it  is  easy  to  show  that- 

«<•'> 

T(-l/r),  =  -y,  t  =  0,l,..., 

T* 


n. 


(3.6) 
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Therefore 

1  &+2) 

2  =  -  ,  i  =  n  -  2,  n  -  3, . . . ,  0.  (3.7) 

*«+i 

The  normal  conditions  in  (3.5)  may  now  be  expressed  as 

4++i1)  #0,  i  =  n-2,n-3,...,0.  (3.8) 

Analogous  to  [6],  [11],  and  [19],  we  then  have 

THEOREM  3.1.  The  polynomial  F(c)  6  S[c]„  is  stable  iff 

I.  ^  0,  t  =  n  —  1,  n  —  2, . . . ,  1,  and 

II.  «/„  =  vzr{T(0)n,T(0)n_u...,T(0)o}  =  0. 

One  of  the  main  advantages  of  the  complex-#- BT  is  that  all  computations  may  be 
carried  out  through  real  arithmetic  only  [19].  The  same  holds  true  for  the  the  complex-6- 
BT  introduced  above  as  well.  To  see  this,  let 

T(c)i  =  S(c)i  +  jA(c)i  with  6,  =  Re[6,]  +  jlm[<5,],  (3.9) 

for  %  =  2, 3, . . . ,«.  It  is  easy  to  show  that  5(c),- ’s  and  A(c);’s  form  sequences  of  6-symmetric 
and  6-antisymmetric  polynomials,  respectively.  Now,  (3.3)  may  be  expressed  as 

5(c), _2  =  —  [Re[6,](2  -I-  rc)  •  5(c),-!  +  Im[<5;]rc  •  A(c);_i  -  5(c),] ; 

1  +  rc  (3  10) 

A(c),_2  =  — - [— Im[6,]rc  •  5(c),_i  +  (2  +  rc)Re[6;]  ■  A(c),_i  -  A(c),-] , 

1  +  rc 

for  i  =  2, 3, . . . ,  n. 

Remark.  Note  that,  T(0);  =  5(0);  +  j'A(0),  =  5(0);. 

In  the  real-6-BT  construction,  a  certain  ‘scaling’  of  {T(c), }•*_(,  was  useful  [6],  We 
use  the  same  technique  in  the  complex-6-BT  case  as  well,  thus  providing  the  following 
advantages:  (a)  Terms  containing  r  are  avoided  during  construction,  (b)  6,  and  »/,-  may  be 
deduced  by  simple  inspection,  and  thus  (c)  computational  effort  is  reduced. 
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The  sequence  of  polynomials  that  incorporates  ‘scaling’  is  { C/ ( C ) * } ”=  o  where 

=  ,,,<=> «*’  =  (-7)  4'’. <■•  =  0,1 — «, 

Jfc=0  c  </r  \  / 

for  *  =  0, 1, . . .  ,n.  Thus,  from  (3.3),  we  get,  for  i  —  n  -  2,n  —  3, . . .  0, 


(t)  /  £  .  7  \  (i+l)  («+2) 

u0  =  (^i+2  +  ^*+2)^0  —  W0  i 


«i°  =  (<i+2 + <i+2)u<l'+1)  -  - «S7« + «i" „  *  =  1,2,. 

Note  that 

1  1  1 

2- 3- 

T  *i+i  ui+ 1 


.(*'+»)  ..(i+2)  ,  ..(<) 


(3-11) 


(3.12) 


(3.13) 


and 

un  =  var{T(0),}*=o  =  var{uj°}f=0.  (3.14) 

Therefore,  condition  II  of  Theorem  3.1  may  be  checked  by  inspecting  the  constant  coeffi¬ 
cients  of  { 17(C).  }”=o- 


Remark.  One  may  use  the  same  ‘scaling’  strategy  in  an  implementation  that  uses  only  real 
arithmetic. 


Relationship  between  complex-q-BT  and  complex-8 -BT 

As  was  agreed  upon  previously,  given  .F(c)n  e  $[z],  let  us  use  the  notation  F(z)n  to 
indicate 

F(z)„=\F(c)n  I  (3.15) 

I  C — *  Z 

where  A  G  5R  is  a  possible  scaling  constant.  The  establishment  of  the  relationship  be¬ 
tween  the  rows  of  complex-^-BT  of  F(z),  i.e.,  {T(z),}"_0,  and  complex-6-BT  of  F(c ),  i.e., 
(T(c),}"_0,  which  is  the  subject  of  this  section,  is  useful  later  in  obtaining  the  Schur-Cohn 
minors  from  the  latter. 

Cl-AIM  3.2. 


Ft(z)n  =  AF*(c) 


c.  — *  z 
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Proof.  Note  that 

F'(z)„  =  z-'P  (I)  =  \z"F(c)n  |^=  \z"F  (L^±)  ; 

f*(c)„  =(1+’'c)“;'(-nf-)  I  =  — )  ■ 

c^z  \  1  +  TC/n  lc— 2  \  TZ  J  n 

The  claim  is  thus  proven.  ■ 

Theorem  3.3.  The  rows  of  the  complex-?- BT  of  F(z)  and  the  complex-<$-BT  of  F(c)  are 
related  by 

s  I 

,  i  =  n,n  -2,...; 


n*)i  =  < 


AT(c)i 

;T(c), 


,  *  =  n  —  l,n  —  3, _ 


Proof.  First,  using  Claim  3.2,  note  that 

T(z)n  =  A  T(c)n 

Thus,  Theorem  3.3  is  established  for  t  =  n.  i  =  n  —  1  may  also  be  established  directly. 
For  i  =  n  —  2,  n  —  3, 

Corollary  3.4. 


,0,  use 

(2.3)  and  (3.3). 

II 

r 

l  r«  ti  i 

for  i 

i  a  t(,) 

for  i 

II 

f 

I  rt S  i 

for  i 

L  r'  +  l  C»  » 

for  i 

Proof.  This  follows  directly  from  Theorem  3.3. 


4.  Schur-Cohn  Minors  for  6-Systems 


We  now  develop  quantities  that  may  be  considered  the  analogs  of  Schur-Cohn  minors  for 
6-system  polynomials. 

LEMMA  4.1.  The  relationship  between  the  complex-6-BT  of  F(c)n  6  5[c]„  and  the  Schur- 
Cohn  minors  A;,  t  =  1,2, ...  ,n,  of  F(z)n  6  S[z]„  is 

A2 


A,  = 


2  r2( 


\ln— i+1  ln—i  '  ln  —  «+l  ln  —  i  —  \ 


A2 


2  r2( 


] ,  A0  =  l,  A<  =  o,  t  <  o. 


Proof.  Note  that,  the  relationship  between  the  complex-g-BT  of  F(z)n  and  its  Schur-Cohn 
minors  are  given  by  [25] 


A  1 

A*  ~  2 


/if n— n— «)  rfn-«+lWn-»)\  *  1  |Wn-i+l),^n-.)|2X 

tf0  tn-i  +  Si-i+1  t0  )&x~l  ~  ~\tn_i+1  tn_i  |  A,_2 


with  Ao  =  1  and  A,  =  0,  i  <  0.  Now,  the  claim  follows  from  Corollary  3.4. 
Let 

Then,  from  Lemma  4.1,  A  in  Theorem  2.5  is  given  by 


(4.1) 


A  =  A2  •  D  ■  A  ■  D 


(4.2) 


where 


A  = 


Re[t(„”)t{nn_-1)] 

2^n-l  ln  — 2  J 

0 

0 

r  r#( «  —  1 )  .-(n-2)i 

2  lln  — 1  n-2  J 

0 

0 

T  f/(n-2)j('*-3)l 

2  —  2  ln-3  J 

0 

(4.3) 

0 

0 

0 

Rc[4,)C)1- 

Clearly,  positive  definiteness  of  A  and  A  are  equivalent  statements.  Hence,  we  may 
consider  the  principal  minors  of  A  to  be  the  Schur-Cohn  minors  of  F(c). 
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DEFINITION  4.1.  The  Schur-Cohn  minors  of  F(c)  G  S[c]„  are  the  principal  minors  of  the 
tridiagonal  Hermitian  matrix  A  in  (4.3). 

Therefore,  from  Theorem  2.4,  we  have 

THEOREM  4.2.  The  polynomial  F(c)  €  S(c]n  is  stable  iff  A,  >0,  i  =  1,2, .. .  ,n,  where  A j 
is  the  (t  x  t)-principal  minor  of  A  in  (4.3). 

Remarks. 

1.  Tridiagonal  Hermitian  matrices  constitute  an  important  class  of  matrices  that  have 
been  extensively  investigated  in  matrix  theory  literature  [26].  See  also  [10]. 

2.  Since  the  Schur-Cohn  minor  A i  obtained  from  the  complex-g-BT  are  necessarily 
proper  [10],  [25],  the  Schur-Cohn  minors  defined  above  for  6-systems  are  proper  as  well. 

In  terms  of  the  ‘scaled’  sequence  of  polynomials  {!/(C)i}”=o>  Theorem  4.2  may  be 
stated  as 

COROLLARY  4.3.  The  polynomial  F(c)  e  S[c]„  is  stable  iff  A,  >  0,  t  =  1,2,  ...,n, 
where  A,-  is  the  (i  x  t)-principal  minor  of 


2  tun  —  1  un- 2  i 

o 

0 

2  i  n  —  1  Un-2  J 

-ReK_,  K-2  1 

••• 

0 

0 

2  —  2  Wn- 3  J 

0 

(4.3) 

0  0  0  •  •  ■  — Re[ujl^Up°^] 


Proof.  Using  (3.11),  and  factoring  out  the  appropriate  diagonal  matrices,  the  result  im¬ 
mediately  follows.  ■ 

Remark.  Again,  notice  how  the  use  of  the  ‘scaled’  sequence  simplifies  the  entries. 


5.  Algorithm  for  Checking  Stability  of  2-D  6-Systems 

To  check  condition  II  of  Theorem  2.2,  we  may  adopt  the  following  approach: 

(a)  Express  F(ci,c2)  €  ^(ci]n, [c2]nj  as  a  polynomial  in  S[c2]nj  so  that  its  coefficients,  as 
well  as  the  corresponding  Schur-Cohn  minors,  are  parameterized  by  ci  G  T$.  Here,  we 
have  assumed  that  ni  >  n2;  otherwise,  the  roles  of  ni  and  n2  may  be  interchange* 

(b)  Check  positivity  of  each  of  the  Schur-Cohn  minors,  or  positive  definiteness  *  „ne 
tridiagonal  Hermitian  matrix  A  6  Sn,Kni,  for  all  ci  6  T<  (see  condition  II  of  Theorem  2.2 
and  Theorem  4.2).  These  checks  may  be  simplified  by  applying  a  direct  extension  of  Siljak’s 
result  [16], 

However,  construction  of  the  complex-6-BT  and  the  entries  of  A  require  complex 
conjugation  of  certain  entries  that  are  functions  of  cy  6  T&.  This  of  course  complicates  the 
scheme  since  cj  =  — Ci/(1  4-  rc i),  Vex  €  Ts-  On  the  other  hand,  in  dealing  with  2-D  q- 
system  stability,  we  have  zj  =  1/zj,  Vzi  €  Tq.  This  simple  relationship  has  led  to  stability 
checking  schemes  that  use  the  complex  forms  of  tabular  forms  [10]  that  incorporate  the 
polynomial  array  method  [27].  To  circumvent  the  above  difficulty,  the  algorithm  given 
below  uses  the  real-6-BT  in  order  to  check  Theorem  2.3.  In  the  appendix,  an  easily 
implementable  algorithm  that  yields 

G(x,c2)  =  G(ar)ni(c2)2nj  =  F(c1,c2)F(c1,c2)|  c,6r,  G  3?[x]„l[c2]2„2  (5.1) 

x=(«l+«l)/2 

is  provided.  Note  that 

ci  €  x  G  [— 2/r,0].  (5.2) 

Before  proceeding,  however,  it  is  important  to  note  that  tabular  methods  are  useful 
in  checking  for  no  roots  to  be  outside  the  stability  region.  However,  since  in  typical  2-D 
stability  studies  the  2-D  transforms  are  taken  with  positive  powers  [9-10],  prior  to  applying 
the  stability  check,  the  following  ‘preparation’  must  be  done: 

(a)  Condition  I  in  Theorem  2.3  may  be  checked  by  explicitly  finding  the  roots  or  applying 
the  real-6-BT  to  ensure 

F8(ci)(-lA)  =  (l  +  rc1)"'F^r^-)  (-^j  /  0,  Vc,  (5.3) 
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(that  is,  polynomial  is  reciprocated  with  respect  to  Ci). 

(b)  First  form 


G(*)„l(ci)j„,  =  ^TV,2ni)(x)c2  G  *[*!«» [ca]j«,  where  g^Jni)(z)  =  €  »[*]»,.  (5.4) 

1= 0  fc=0 

Here  x  €  (— 2/r,  Oj.  Now,  condition  II  in  Theorem  2.3  may  be  checked  by  applying  the 
real-6-BT  to  ensure 


2  nj  nj 

G(*)n,(c2)2n,  =  ^S(/nj)(*)c2  where  j[2n,)(x)  =  €  »[*]«, 

t=0  i=0 

=  G(*)'(C2) 

=  (l+rc2)2niG(z)(r^-)  *0,  Vx  6  [— 2/r,  0],  Vc2  €  9  \  Us 


(5.5) 


(that  is,  polynomial  is  reciprocated  with  respect  to  C2).  Again,  x  €  [— 2/r,0|. 


We  will  hence  implicitly  assume  that  the  given  2-D  6-polynomial  has  already  been  ap¬ 
propriately  ‘prepared’  as  above.  In  addition,  the  construction  of  the  real-6-BT  for  G(x)(c2) 
requires  ensuring  (11] 

9o2n2\x)  £  0  and  9?^  >  ^x  6  [~2/r,0].  (5.6) 

Violation  of  the  first  condition  in  (5.6)  is  equivalent  to 

F(ci)(0)  =  0  for  some  cx  6  T«.  (5.7) 

Assuming,  with  no  loss  of  generality,  >  0  for  some  x  €  (— 2/r,  0],  violation  of  the 

second  condition  in  (5.6)  is  equivalent  to 

F(cx)(  —  1/r)  =  0  for  some  cx  €  T&.  (5.8) 

Therefore,  each  of  these  violations  imply  instability.  Verifying  condition  (5.7)  must  be 
included  in  the  algorithm.  Condition  (5.8)  is  automatically  verified  when  condition  I  in 
Theorem  2.3  is  checked  (see  (5.3)). 
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Then,  we  have  the  following 

Theorem  5.1.  The  2-D  6-system  in  (2.7)  is  stable  iff 

I.  F(ci)(-l/r)  ^  0,  Vci  G  Us,  and 

II.  F(ci)(0)  /0,  Vc,  G  Ts,  and 

III.  A,(0)  >0,  Vi  =  1,2, . . .  ,2n2,  and 

IV.  A2n,(a:)  >  0,  Vx  G  [— 2/r,0),  which  is  satisfied  whenever  A2rti(x)  /  0,  Vx  G  [— 2/r,  0), 
together  with  condition  III. 

Here,  A  is  the  Hermitian  matrix  mentioned  in  Theorem  4.2  corresponding  to  G(x)(c2) 
where  x  G  [— 2/r,  0]. 

Conditions  I  and  II  in  Theorem  5.1  are  easy  to  carry  out  (they  may  in  fact  be  verified 
by  explicitly  finding  the  roots).  Condition  III  and  IV  require  construction  of  the  real-6-BT 
and  the  Schur-Cohn  minors  for  which  we  now  develop  polynomial  arrays  [27].  We  also 
provide  a  scaling  scheme  so  that  the  numerical  reliability  of  the  resulting  algorithm  is 
enhanced. 


5.1.  Polynomial  array  for  entries  of  real-6 -BT 
Express  G(x)(c2)  as 

C(x)(c2)  =  x<n')T  .G.<42"l)  (5.9) 

where  x("‘>  =  [x"‘ ,  x"*"1, . . . ,  1]T,  c<2n’}  =  [c2"’,^"’'1, . . . ,  l]r,  and  G  =  {</,,,}  G 
^(rn+i)x(2n2+i)  jg  coefficient  matrix.  Then,  it  is  easy  to  show  that  [6] 

G(x)(c2)  =  x(n,)T  •  G  •  cfnj)  where  G  =  Gr(2n,)~lp(2niM2n*>.  (5.10) 


Here 

t(2",)  =s  diag{r2nj,r2nj_I,...,l}  G  ft(2n2+i)x(2n2+i). 

P(2ni)  =  {p0}  G  R(2«»+1)x(2u1  +  I)  where  Pi.  =  +  l -ipi. 

The  elements  pij ,  which  in  fact  are  those  of  the  Pascal’s  triangle,  are  given  by 

(  0,  for  i  <  j; 

1,  for  i  —  j] 

Pi- i,j-i  +  Pi-i,),  elsewhere. 


(5.11) 


Pxj  = 


(5.12) 


•  I 
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The  real-5-BT  is  constructed  using  the  ‘scaled’  polynomial  sequence  in  (3.11-14).  Let 


H(y)(0  =  G(x)(c2)  | 

HWO  =  G(z)‘(c2)  |  =  G(x)(ct)  |  . 

*  =  -*/«•  *  =  -»/r 


(5.13) 


Note  that,  x  G  [-2/r,0j  iff  y  6  [0,2].  Now,  using  (5.9-12),  row  #2n2  and  2n2  -  1  of  the 
corresponding  ‘scaled’  real-6-BT  are  given  by 

2n, 

e~o 

—  y(ni)T  .  f(ni)_lQT(2nl)-l^j(2ni)  +  p(2n,)^  _  ^{2n,). 

u(y)((h*,-i  =  e‘  -  ^)(C)7/gfa)(C)  (514) 

U  ~UT 

—  y(”‘)T  .  •j^n‘)-lrG7J(2n*)"l(/(2n,)  _  p(2n2)^  ^ 

where  C*2n,)  =  [C2nj,  <2ni-\ . . . ,  1]T,  and 

T*ni)  =diag{(-r)n*,(-r)n‘“1,...,l}  G  ft(n>+1)*(n‘+1>; 

J(2nj)  =  diag{(— l)2n*,  (—l)2”2-1, . . . ,  1}  G  £(2n’+1)x(2n2+1>;  (5.15) 

P(2nj)  =  {P0}G*(2ns+1)x(2n,+1)  where  Pij  =  {-l)i+j  Pij. 

Each  element  of  the  remaining  rows  is  of  the  form 

( «)/  s 

«i')(y)=  <  =  0,1,...,*,  t  =  2n2,2n2  -  1,...,0,  (5.16) 


(5.15) 


(5.16) 


where  n(,°(y)  G  &[y]„(o  and  d(,)(y)  G  9R[y]c<.).  Substituting  in  (3.12),  it  is  easy  to  show 
that,  for  t  =  0, 1, . . .  ,£, 


-4”  =  ■SS’W  -  2n^+ ")  -  nK"n}‘«>  +  n«,  for  i  =  2„,  -  2,. . .  ,0; 


<0  “  {  rff'+!)n' 


for  i  =  2n2, 2n2  —  1, 
1+4,0’  for  i  =  2n2  —  2, . . , ,  0. 


(5.17) 


Note  that  u^2"2)  =  n^2'1^  and  u[2n7  =  7i^2n’  Moreover 


for  t  =2n2,2n2  -  1, 

(cr<,+2)  4-a<,+1),  for  i  =  2n2  -  2, . . . , 0; 

(.)  _  |  0,  for  i  =  2n2,2n2  -  1, 

(  cr(')  —  n,,  for  i  =  2nz  -  2, . . . ,  I 


(5.  IS) 
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(5.19) 


Scaling  scheme.  Let  us  scale  rows  #2rt2  and  #(2ri2  —  1)  so  that  each  coefficient  takes 
values  in  (—1, 1).  Correspondingly,  for  £  =  0, 1,. . .  ,i;  t  =  2«2,2n2  —  1,  let 

n<‘>  = 

d<*>  =  7<‘Vi>, 

where  A^.t^  >  0,  t  =  2n2,2n2  —  1,  are  the  scaling  constants  and  [:]  denote  scaled 
quantities.  Note  that 

n(2n,)  A(2„,)  *(2n,) 

d( 2n*>  ~  J(2nj)  ’ 

rc(2na-l)  _  A(2n,-1) 

d( 2nj  — 1)  <y(2n2  — I)  ^(2ni-l)  ' 

Now,  substituting  in  (5.17-18),  we  get 


(5.20) 


n 


(2nj— 2) 


A(2n2)A(2n2-l) 
d(2n7-2) 


_ -(2n2),  (2n2-l)  9  -(2n2-lK  (2n2-l)  (2n2) 

-X)  “  n2n,  tn<-l  -  Znl  )  ~  n2nj-l  nt  + 


n 


(2nj  —2) 
/-I 


A(2’**)A(2n*-1>  ’ 


(5.21) 


_  -a  n2  i  . 


(5.22) 


7(2nj)A(2n2-l)  "  2n*' 

It  can  now  be  seen  that,  it  is  only  necessary  to  compute  the  quantities  on  the  left  hand 
side  of  (5.21).  Then,  one  may  scale  these  to  get 

n(2n2-2)  _  A(2n2)A(2n2-I)A(2n2-2)^(2n2-2). 

^2n2-2)  _  7(2n2)7(2n2-2)A(2n2-I)J(2n2-2) 

Note  that 

n(2n2-2)  A(2nj)A(2,,J-2)  ri*2"2-2) 

d( 2n2-2)  =  7(2n2)7(2nj-2)  J(2n7-2)  '  (5.23) 

Continuing  in  this  manner,  the  computation  of  the  entries  of  real-6-BT  may  be  summarized 
as  follows: 

(a)  From  (5.14),  compute  n^\d^\  i  =  2ri2,2ri2  —  1. 

(b)  From  (5.19),  use  scaling  constants  A^,  7^,  i  =  2ri2,2ri2  —  1,  to  get  i  = 

2n2, 2n2  —  1. 

(c)  From  (5.21),  for  t  —  0, 1, ...  ,  t;  i  —  2ri2  —  2,. . .  ,0,  compute 


n 


(0 


(,)  —  ni+ 2  K-1  ^ne  )  n«+i  nt 


(_  _  *(«+2)^<.+l) 

K(n 

^(,)  _  Jb+2) -(«+*) 
tA>)  ~a  n«+>  ’ 


(<) 

s(«+l)x  ^(«  +  1)a(«+2)  n(-J 


id"  ’ 


(5.24) 


I 
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and  use  scaling  constants  t  =  2n2  —  2, ...  ,0,  to  get  t  =  2n2  —  2, . 

Here,  Kn^  and  axe  constants. 

(d)  Notice  the  relationships 


,0. 


n 


(0 

l  _ 


d(0 


•y(*»*)7(*»i-*)...7(0  3(iT  ’  fort  —  2ri2,2ri2  2, 

A(»»*-1)A<*.,-»)„.A<0  *(0 

tfj) >  tor  t  =  2n2  -  l,2n2  -  3,. . . ,  1. 


5.2.  Polynomial  array  for  Schur-Cohn  minors 

Each  Schur-Cohn  minor  obtained  from  the  table,  in  general,  will  be  of  the  form 

A<(y)  =  6  *(y),  t  =  1, 2, . . . ,  2n2, 

where  N^\y)  €  3f^[y]p<»)  and  £>(,)(y)  €  9fc[y]c<i)-  From  Corollary  4.3,  we  get 


(5.25) 


(5.26) 


A  .{i,\  ,.(2»*2-*+l)  (2na-i)  *  *  unj-'tU'  .l*Bt-'r  A  ;  _  -i  o  / c  OTA 

A*(y)  -  u2m-i+l  U2nj-i  A«-l  Au2n2-t+l  U2  n2-i  A,_2,  *  —  1,  2,  .  .  .  ,  2n2,  (5.27) 


1..(2n*-«+l)a..(2n,-t)* 


where  A0  =  1  and  A ;  =  0,  Vi  <  0. 

Remark.  Actually,  as  in  [10],  one  may  show  that,  for  stability  determination  purposes,  only 
the  numerator  polynomials  of  A;  need  be  computed.  However,  to  contain  the  orders  of  the 
resulting  polynomials,  and  hence  improve  numerically  conditioning,  we  do  not  recommend 
this  scheme. 


Scaling  scheme.  Due  to  the  scaling  of  entries  of  the  real-6-BT,  computation  of  A,,  i  = 
1,2, ...  ,  2n2,  may  be  modified  as  follows:  Let 


A, 


— u 


(2n2) 

2 


U 


(2n,  — 1) 
2n2-l 


-y(n22),y(2n2  — 1)  ^(2nj)^(2n2  —  l) 


(5.28) 


Hence,  it  is  only  necessary  to  compute 


A,  =  — 


-  (2"2) ,-(2n2 
n2n2  u2n2 
J(2n2)  <f(2n* 


(5.29) 
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Continuing  in  this  manner,  the  computation  of  the  Schur-Cohn  minors  may  be  summarized 
as  follows:  From  (5.27),  for  i  —  1, 2, ... ,  2n2,  compute 


A,  = 


.(2n2-i+l)  .(2nj-i) 
n2n2-i+l  n2ni~i 
J(2n2-i+l)J(2n2-«) 


A,_i  +  ~ 


1  \(2nj-«)  -(2nJ-,+  l)^(2n}-,) 

1  n2na-i+l  n2  n2-»  X 

4 -y(2 «J-«)  J(2nj-i+l)  J(2nj-i) 


A, -2 


(5.30) 


where  Ao  =  1  and  A;  =  0,  Vt  <  0. 

Remark.  Note  that,  since  A i(y)  is  necessarily  a  proper  polynomial  (that  is,. denominator  di¬ 
vides  numerator  properly  with  no  remainder),  and  not  a  rational  polynomial  (see  Remark  2 
after  Theorem  4.2),  it  is  easy  to  see  that  S2nt~t+1^S2n,~^  must  divide  n2n" -"i+j1  ^2n" 
exactly. 

5.8.  Algorithm 


The  following  result,  which  is  the  basis  of  the  stability  checking  algorithm,  is  now  obvious 
from  [10]  and  Theorem  5.1: 

THEOREM  5.4.  The  2-D  5-system  in  (2.7)  is  stable  iff 

I.  F(clt  -1/t)  ^  0,  Vex  €  Us,  and 

II.  F(ci)(0)  7^  0,  Vci  €  Ts,  and 

III.  A{(0)  >0,  V*  =  1,2,...,  2n2,  and 

IV.  A2nj(y)  ^  0,  Vy  €  [0,2]. 

The  2-D  stability  checking  algorithm  may  now  be  summarized  as  follows: 

Given. 

A  2-D  5-polynomial  F(ci,c2)  €  3?[ci]„, [c2]nj.  Without  any  loss  of  generality,  assume 
that  ni  >  n2,  and  express  F(cj,c2)  as  F(ci)n,(c2)n,. 

STEP  I.  Condition  I  of  Theorem  5.4: 

Apply  an  explicit  root  location  procedure.  If  result  is  satisfactory,  proceed;  otherwise, 
system  is  unstable. 


STEP  II.  Condition  II  of  Theorem  5.4: 
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Apply  an  explicit  root  location  procedure.  If  result  is  satisfactory,  proceed;  otherwise, 
system  is  unstable. 

Step  III. 

Form  C?(y)(c2)  using  the  algorithm  in  the  appendix;  then  form  U(y)( £)2n}  and 
tf(v)(0*»,-i  from  (5.14).  These  yield  and  Of  course,  d*2nj*  =  d(2n,~1*  = 

1. 

From  (5.19),  obtain  n,2nj\  n^2"J_1\  and  the  associated  scaling  constants  A^2nj) 
and  A*2"5-1).  Of  course,  S2n2^  =  S2n7~^  —  1  and  7(2nj)  =  ■y(2ni~1)  —  i. 


STEP  IV.  Condition  III  of  Theorem  5.4: 

Form  Ai(y)  from  (5.30)  and  check  whether  At(0)  >  0. 

If  result  is  satisfactory,  form  n^2"J~2^  and  S2n7~2^  and  the  associated  scaling  constants 
-\(2nj-2)  and  j(2n*~2)  from  (5.24).  Form  A2(y)  from  (5.30)  and  check  whether  A2(0)  >  0. 
If  result  is  satisfactory,  proceed  likewise  until  A2n2(0)  >  0  is  checked.  Note  that,  this  re¬ 
quires  checking  of  only  the  constant  coefficients.  If  result  is  satisfactory,  proceed;  otherwise, 
if  the  check  fails  at  any  i  =  1,2,...  ,2n2,  system  is  unstable. 

STEP  V.  Condition  IV  of  Theorem  5.4: 

Apply  an  explit  root  location  procedure  to  check  whether  A2nj(y)  /  0,  Vy  6  [0,2). 

Remarks.  The  possible  numerical  difficulties  that  may  arise  in  using  explit  root  location 
procedures  may  be  avoided  as  follows:  (a)  Steps  I  and  II  may  be  verified  using  the  real-6- 
BT  [6],  and  (b)  step  V  may  be  verified  by  the  Sturm  sequence  method. 


6.  Example 


The  stability  checking  algorithm  presented  in  the  previous  section  is  now  illustrated 
through  an  example.  Polynomial  entries  are  denoted  using  a  self-explanatory  shorthand 
notation  where  the  highest  degree  coefficient  is  written  first.  Moreover,  only  four  decimal 
digital  on  the  mantissa  are  shown. 

Consider  the  2-D  polynomial 

'1  50  740  ]  \ 4' 

F(cuc2)  =  [c?  cx  1]  52  2700  38480  c2 

[  740  38480  547600  J  [  1 

with  the  sampling  time  r  =  0.1  s. 

STEP  I.  Condition  I  of  Theorem  5.4: 

By  applying  an  explicit  root  location  procedure,  one  can  show  that 

F(c\ )( — l/r )  =  340c?  +  16680Cl  +  236800  ^  0,  VCl  e  U6. 

STEP  II.  Condition  II  of  Theorem  5.4: 

By  applying  an  explicit  root  location  procedure,  one  can  show  that 

F(ci)(0)  =  740c?  +  38480ci  +  547600  ^  0,  VCl  G  T6. 

STEP  III.  Using  the  algorithm  in  Appendix,  we  get 

1.2800e+03  1.2992e  +  05  5.1904e  +  06  9.6141e  +  07  7.0093e  +  081  r<2' 

G(y)(0  =  (y2  y  1]  5.2480e  +  04  5.4011e  +  06  2.1662e  +  08  3.9968e  +  09  2.8738e  +  10  < 

5.4760e  +  05  5.6950e  +  07  2.2912e  +  09  4.2143e  +  10  2.9987e+  11  J  1  . 

After  scaling,  rows  #4  and  #3  are  computed  as  follows: 

nj4)  =  [1.2859e  -  02,  -5.0693e  -  02,  5.1315e  -  0°’ 
n^4)  =  [— 7.9833e  -  02,  3.1991e  -  01,  -3.2798e 
n^4)  =  [1.9671c  -  01,  -7.9909e  -  01,  8.2798e  -  01); 
n\i]  =  [— 2.3375e  -  01,  9.5836e  -  01,  -l.OOOOe  +  00); 
n(04)  =  [1.1687e  -  01,  -4.7918e  -  01,  5.0000e  -  01], 
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with  A<4>  =  1.1995e  +  12,  and 

n(3J)  =  [— 2.4050e  -  02,  9.4053e  -  02,  -9.4595e  -  02]; 
n^3)  =  [1.3044e  -  01,  -5.1542e  -  01,  5.2252e  -  01]; 
n(!3)  =  [— 2.4703e  -  01,  9.8194e  -  01,  -l.OOOOe  +  00]; 

43>  =  (l-6469e  -  01,  -6.5463e  -  01,  6.6667e  -  01], 
with  A^3)  =  5.3490e  +  10.  Of  course,  d^4)  —  <?3)  =  1  with  7I4)  =  7O)  =  1. 

STEP  IV.  Condition  III  of  Theorem  5.4: 

We  get 

Ax  =  [3.0926e  -  04,  -2.4286e  -  03,  7.2184e  -  03,  -9.6217e  -  03,  4.8541e  -  03], 

Clearly,  A^O)  =  4.8541e  -  03  >  0. 

Now,  row  #2  is  computed  as  follows: 

n^2)  =  (— 8.8619e  -  03,  6.8253e  -  02,  -1.9920e  -  01,  2.6099e  -  01,  -1.2957e  -  01]; 

n(x2)  =  [3.3584e  -  02,  -2.5969e  -  01,  7.6068e  -  01,  -l.OOOOe  +  00,  4.9793e  -  01]; 

*o2>  =  [-3-3584e  -  02,  2.5969e  -  01,  -7.6068e  -  01,  l.OOOOe  00,  -4.9793e  -  01], 

with  A(2)  =  4.2420e  -  02.  Also, 

d{2)  =  [— 2.5424e  -  01,  9.9428e  -  01,  -l.OOOOe  +  00], 
with  7^)  =  9.4595e  -  02.  We  get 

A2  =  [1.8046e  -  07,  -2.8190e  -  06,  1.9343e  -  05,  -7.6148e  -  05,  1.8810e  -  04, 

-  2.9857e  -  04,  2.9737e  -  04,  -1.6992e  -  04,  4.2654e  -  05], 

Clearly,  A2(0)  =  4.2654e  -  05  >  0. 

Now,  row  #1  is  computed  as  follows: 

*1°  =  (2  5168e  -  03,  -2.8515e  -  02,  1.3555e  -  01,  -3.4597e  -  01,  5.0000e  -  01, 

-  3.8792e  -  01,  1.2623e  -  01]; 

”oU  =  {— 5.0336e  -  03,  5.7031e  -  02,  -2.7110e  -  01,  6.9194e  -  01,  -l.OOOOe  -f  00, 
7.7584e  —  01,  -2.5246e  -  01], 
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with  A*1*  =  3.0980e  -  02.  Also, 

d(1)  =  [— 3.3954e  -  02,  2.6151e  -  01,  -7.6322e  -  01,  l.OOOOe  +  00,  -4.9646e  -  01], 
with  7^  =  2.6099e  -  01.  We  get 

A3  =  (4.0500e  -  10,  -9.3525e  -  09,  9.9260e  -  08,  -6.4020e  -  07,  2.7947e  -  06, 

-  8.6990e  -  06,  1.9797e  -  05,  -3.3188e  -  05,  4.0679e  -  05,  -3.5552e  -  05, 
2.1029e  -  05,  -7.5594e  -  06,  1.2489e  -  06]. 

Clearly,  A3(0)  =  1.2489e  -  06  >  0. 

Now,  row  $0  is  computed  as  follows: 

h<0)  =  [— 1.0487e  -  04,  1.9379e  -  03,  -1.6174e  -  02,  8.0291e  -  02,  -2.6251e  -  01, 
5.9070e  -  01,  -9.2642e  -  01,  l.OOOOe  -j  00,  -7  1104e  -  01,  3.0076e  -  01, 
-  5.7473e  -  02], 

with  A<°>  =  4.4719e  -  02.  Also, 

of*0)  =  [— 6.7946e  -  04,  1.0355e  -  02,  -6.9373e  -  02,  2.6679e  -  01,  -6.4420e  -  01, 
l.OOOOe  +  00,  — 9.7458e  -  01,  5.4519e  -  01,  -1.3404e  -  01], 
with  7^  =  9.4174e  —  01.  We  get 

A4  =  (4.3531e  -  12,  -1.3058e  -  10,  1.8400e  -  09,  -1.6166e  -  08,  9.9118e  -  08, 

-  4.4970e  -  07,  1.5618e  -  06,  -4.2352e  -  06,  9.0628e  -  06,  -1.5355e  -  05, 
2.0530e  -  05,  -2.1433e  -  05,  1.7129e  -  05,  -1.0130e  -  05,  4.1814e  -  06, 

-  1.0762e  -06,  1.3014e-07], 

Clearly,  A4(0)  =  1.3014e  -  07  >  0. 

STEP  V.  Condition  IV  of  Theorem  5.4: 

By  applying  an  explicit  root  location  procedure,  one  can  show  that 

A4(y)^o,  Vye[0,2]. 


Thus,  we  conclude  that  F(cj,C2)  is  stable. 


«  * 


t 


I 


7.  Conclusion  and  Final  Remarks 

In  this  paper,  we  have  developed  an  efficient  stability  checking  algorithm  applicable  for  2- 
D  6-system  characteristic  polynomials.  Our  purpose  here  is  to  obtain  a  direct  algorithm 
due  to  the  possible  numerical  disadvantages  associated  with  indirect  methods  that  utilize 
transformation  techniques. 

In  arriving  at  the  algorithm,  the  following  contributions  have  been  made:  (a)  Tab¬ 
ular  method  of  stability  checking  applicable  for  6-system  polynomials  possibly  possessing 
complex-valued  coefficients,  (b)  quantities  that  may  be  regarded  as  the  Schur-Cohn  minors 
applicable  for  such  systems,  and  (c)  polynomial  arrays  for  computing  both  table  entries 
and  Schur-Cohn  minors. 

The  proposed  Schur-Cohn  minors  lets  one  use  a  Siljak-like  simplification  [16]  in  the 
stability  check.  Although  the  algorithm  utilizes  only  the  real-6-BT,  results  regarding  the 
Schur-Cohn  minors  are  in  fact  valid  for  the  more  general  complex-valued  coefficient  case 
as  well. 

As  in  [10],  it  is  possible  to  develop  the  algorithm  such  that  only  the  numerator  poly¬ 
nomials  of  the  entries  of  the  real-6-BT  and  the  Schur-Cohn  minors  are  computed.  Then, 
we  do  not  require  polynomial  division  operations.  However,  our  experience  has  been  that 
such  a  scheme  is  prone  to  be  numerically  unreliable.  This  is  mainly  due  to  the  explosion 
of  polynomial  degree  especially  in  computing  the  Schur-Cohn  minors.  To  avoid  these  dif¬ 
ficulties  and  enhance  numerical  reliability,  we  have  (a)  introduced  a  scaling  scheme,  and 
(b)  used  polynomial  division  to  contain  the  polynomial  degree.  The  latter  is  not  new;  in 
fact,  MJT  also  uses  this.  If  the  user  is  interested  in  implementing  the  algorithm  using 
PRO-MATLAB  [28],  these  polynomial  division  operations  may  be  conveniently  performed 
using  the  routine  deconv. 

We  believe  that  a  suitable  scaling  strategy  can  improve  the  numerical  reliability  of 
the  MJT  as  well.  The  authors  are  currently  looking  into  this. 

The  algorithm  developed  is  easily  implementable  on  a  computer.  The  authors  have 
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implemented  it  via  a  C-language  routine  that  the  interested  reader  may  request  from  the 
second  author. 
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Appendix.  Algorithm  to  obtain  G(x)ni  (02)20,  from  F(ci)ni(c2)n, 


Given 


2  1 

^(cl)n1(c2)n,  =  Y1  ' 4  where  -MC1 )  =  X  ^k,l  '  Cl»  C»  6  (a*1) 


we  now  develop  an  algorithm  that  yields 


*n2 

G(x)nt  (02)2,1,  ~  ^  ]9j(x)  '  4  =  -^'(f'l)ni(c2)n,  ‘  ■^'(ci)ni(c2)n,  >  C1  6 


First,  we  see 


G(x)(c2)  =  ££;/,(c1)/j(c1)-c'+> 

'=°  (a.3) 

n,  n,+<  2 n,  j  '  ' 

=  £  £  ■  4  =  Y2Y^ /<(ci)/;-<(cj) •  c2 

/=o  i=<  i=o  r=o 

(quantities  with  negative  subscripts  are  taken  to  be  zero).  Hence,  comparing  (a.2-3),  we 
get 

J  j  fn,  n,  1 


9j(x)  =  X  Mc  1  )/y-<(ci )  =  X  XX Mij-t  ■  cid\ 


(=0  Lfc=0  t=0 


where 


}  ni 


=  X  X  f^tfkj-e  •  (CiCi)*  +  X  X  ■  c\c\ 


e=o  k= 0 


]  f»l 


Jt = 0  i=o 
•  ** 


*=0  fc=0 


*  =  X  X  X  fk,tfi,j-t  •  cicl  +  XI  fk,efi,j-e  •  C$c\ 


(=0  k= 0  •=« 


j  n,  k 


ni  k 


=  X  XX  fk,tfij-t  ■  4c\  +  X  X  Mkj-t  ■  c\  4 


t-0  k= 0  •=» 

i*k 


j  Hi  k 


k= 0  •=» 


—  X  X  X  ( fk,(fi,j-t  ■  ci  +  fi,efkj-e  ■  c,  )  •  (cic,)’ 


1=  0  k= 0  -=o 

L  •** 

«i  A- 


=  £  £(«i«. )‘  £  [/m/w-«  •«?-*  +  hi-'h.'  «{-')  ■ 


(a.5) 
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Let  us  use  the  notation 


4">  =  ^L±£L,  c,  €  7i,  n  =  0, 1, - 


Noting  that,  for  Cj  £  7j, 


it  is  easy  to  shov  that 


Substituting  in  (a.5),  we  get 


ci  = 


1  4-  TCi  ’ 


2  (j) 

CiCi  =  — c\  . 


*  =  E 


J  nl  /  o\  1 

E  EE2^'/--'(t) 

/=o  U=o  i=o  v  7 


Substituting  in  (a. 4),  we  get 


J  ni 


9j(z)  =  Y,  S 


(=0  k= 0 


fk,t{k.j-l  ££)‘  +  E2A,/,,-,  (^)‘  ■ 


Now,  in  order  to  develop  the  algorithm,  we  need  a  recursive  procedure 
pute  Cj  ,  n  =  0, 1, . . ..  To  proceed,  we  note  that 


.(*»)  _  (ci+cx)(c"~  +  cf~  )  -cici(c"  -t-c”  l) 

2 


ci  = 


Let 


where 


t=0 


c1,11  =  X. 


Remark.  Note  that 


=  1. 


Substituting  (a.12)  in  (a.ll),  and  equating  similar  coefficients,  we  get 


(n-l) 


1 

+  -c 


T 


(n-2) 

l,t-I 


) 


,  i  =  0,...  ,n,  n 


2,3,. 


(a.6) 

(a.7) 

(a.8) 

(a.9) 

(alO) 

to  com- 

(a.ll) 

(a.12) 

(a.  13) 

(a.  14) 

(a.  15) 


I 
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For  instance,  Cj"\  n  =  0, 1, ...  ,5,  may  be  conveniently  obtained  from 


r1  i  r 1 

0  1  x 

0  2/r  2  I2 

0  0  6/r  4  x3 

0  0  4/r2  16 /r  8  x4 

.0  0  0  20/ t2  40 /r  16 J  lx5 


